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In  the  Calculation  of  Chealcttl  Traneltion  Probabilities 


A Si^wm^ry  Report 


The  work  at  hydrocarbon  Research,  Inc.  laoder  Contract 
E^^OSR.737  with  the  Office  of  Naval  Research  has  been  devoted  to 
the  development  of  transition  probabilities  for  displaoenent  reac- 
tions in  a form  suitable  for  calculation  on  large  scale  automatic 
conqjuting  machinery <,  The  formulation  of  the  transition  probablli'. 
ties  is  based  i^n  the  quantum  theory  of  chemical  kinetics  developed 
previously  try  the  authors  The  present  work  is  part  of  a combined 

experimental  and  theoretical  study  of  hg-Tdrogen-oxygen  flames  carried 
out  Jointly  by  hydrocarbon  Research,  Inco  and  the  National  Bureau 
of  Standards,  and  the  calculations  based  upon  the  model  described 
here  are  to  be  carried  out  on  SEAC,  the  automatic  computar  of  the 
Bureau  of  Standards-,  Final  coding  of  the  problem  for  SEiAC  is  cur- 
rently being  performed  at  the  Computation  laboratory  of  the  Bureau,.-. 


(1)  S,  Golden  and  &.  N.  Falser,  **The  Quantum  Mechanics  of  Chemical 
Kinetics  of  Homogeneous  Gas  Phase  Reactions",  J,  Chem,  Pbys., 
(17)  1949,  pp.  630-643.-  'I'he  model  used  here  differs  from  the 
earllei-  one  in  several  important  ysya,  in  particular  with  re-- 
gard  to  the  treatment  of  the  soatterii^  and  the  chemical 
latent  ials 
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The  iicasdiate  objective  baa  been  to  provide  a meaua  of 
calculating  transition  probabilities  for  the  reactions 

O2  * S — ■ I OB  •>  0 

^2  0 .11 1 1- 1.  OH  ■4’  B;, 

but.  the  formulation  is  sufficiently  general  so  that  it  may  be  ap- 
plied to  a large  class  of  displacement  reactions  „ 

The  pzdncipal  aim  of  the  present  work  has  been  to  arrive 
at  a method  of  computation  that  (1)  would  not  require  machine  capacity 
in  excess  of  that  of  SBAC,  and  (2)  could  be  carried  out  on  SSAC  in 
minimal  time.  We  have  STiceeeded  in  satisfying  the  first  of  these 
requirements'  Precise  statements  concerning  computing  time  cannot 
be  made  until  the  coding  is  completed,  but  the  representatives  cf 
the  Confutation  laboratory  have  sstiioated  that  a single  transition 
(from  a given  vibration-rotation  state  of  the  reactants  to  a given 
vibration-rotation  state  of  the  products)  can  be  computed  in  under 
five  minutes o On  this  basis,  it  appears  that  40  - 50  hours  of 
confuting  time  will  be  required  to  calculate  a given  reaction  at  a 
given  teraperatuire . 

The  formula  for  computing  a given  transition  probability, 
which  has  been  arrived  at  as  t)w  result  r.  prociouTi  studies,  may 


T = 

Uf  m *= 

K,  L e 


ten;|)erature, 

quantUB  cumbers  corresponding  to  vlbxstlon  for 
Initial  and  final  nolaeales,  respectively; 
quantim  nonbers  ooirrespondlng  to  rotation  for 
initial  and  final  laolacules,  respectively; 
noTBalised  Bermlte  polynonial  of  order  n; 
spherical  Bessel  function  of  order 


El,  E2 


Legendre  polynomial  of  order  s; 

<=  discrete  energy  levels  for  in^.tlal  and  final 
molecule,  respectively  • 


“ / /C I ♦“  ^ * ' ) I 


The  quantities  Bj,  aj^,  bj^,  c^,  1 = 1,  2^ denote  constants  which 

are  determined  from  the  physical  properties  of  the  system..  Defini- 
tions of  these  constants,  together  with  their  mnoerlcal  values  for 
the  hydrogen^oxyfTen  reaction,  are  given  in  Appendix  Ic  The  quanti- 
ties Ag,  s = 0,  1,  2,  denote  the  coefficients  in  the  expansion 
in  a Legendre  series  of  the  pertiirbation  potential  assumed  responsible 
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for  chemical  trazisltlon,  naioelyy 


where  ^ denote  the  vectors  joining  the  atoms  of  the  initial 

and  final  moleooleSy  respectively.  Numerical  values  of  the  A's  are 
given  in  Appendix 

In  the  eases  of  interest,  it  has  been  found  that,  to  a good 
Epprcximation,  one  m^  take  Ag  -=  0 for  s>3.  This  assumption  leads 
to  considerable  simplification  in  the  calculations,  for  than  the 
quantities  vanish  unless  |k  - 1*|^  3,  |l  - 3c- 

FroB  the  e<mputational  point  of  view,  the  two  major  dlf  ^ 
floultlas  are  the  calculation  of  the  double  Integt^l  C,  and  the  evalua- 
tion of  the  appropriate  Bessel  functions.-.  We  shall  consider  these 
in  turn.. 

Various  analytical  methods  vei*e  employed  in  an  effort  to 
reduce  the  double  integral  in  G to  a convenient  computing  form,  uut 
we  were  unable  to  obtain  az^  suitable  results...  It  seems  unlikely 
that  one  could  find  a closed  form  foi  G,  and  the  various  efforts  to 
expand  in  a aeries  led  to  slowly  converging  terns  of  very  great  conn- 
plexlty.  Accordingly,  it  bar  been  decided  to  evaluate  G directly 
by  means  of  nuomerloal  quadrature  c 
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It  will  be  seen  t>iat  this  approach  places  the  burden  of 
the  vorlc  on  the  calculation  of  large  numbera  of  3essel  functions-. 
Thus,  it  becomes  iiiqx>rtant  to  use  a qtiadrature  fomtula  which  requires 
the  fewest  possible  number  of  points.  In  view  of  the  weight  function 

expr*-  f y*'  which  appears  In  the  Integrand,  the 

appropriate  integration  formula  to  use  here  is  the  Gauss-type  quadra- 
ture based  upon  the  Renoite  polynomials.^^) 

The  optimum  number  of  points  to  be  used  in  this  quadrature 
has  been  estimated  by  examinatjon  of  the  errors  obtained  in  the 
■erlcal  integration.  For  integrands  with  parameters  in  the  range  of 
interest  of  the  hydrogen-oxygen  reactions,  and  for  transitions  up 
to  the  third  excited  rlbraticnal  level  (m,  a = 0,  1,  2,  3),  we  have 
found  that  the  error  resulting  from  a six-point  quadrat\u«  formula 
in  each  variable  is  veil  less  than  half  of  one  percent.  A flvo~ 
point  formula  waa  fcimd  to  yield  errors  up  to  ten  percent,  particit- 
larly  at  the  higher  vibrational  levels..  Accordingly,  the  double  in- 
tegrals in  G will  be  evaluated  by  means  of  the  six-point  formula^ 
thus  reqtiirlng  that  each  ii.tegrand  be  evaluated  at  36  points  c 


(2)  H,  E.  Selzer,  R.  Zucker,  and  R.  Capuano,  **Tabla8  of  the  Zeros 
and  Wel^t  Factors  of  the  First  Twenty  Hermlte  Polynomials'*, 
J.  Research  Kat.  Bur.  Stds..,  (48)  1952  pp.  111-116. 


Tho  final  Intogratiooi  for  the  transition  probability  will 
be  performed  by  means  of  a Gauss^type  quadrature  based  upon  the 
laguarre  polynomials  After  euitable  rearrangement,  the  transi- 

tion probability  will  then  be  calculated  on  SEAC  in  the  following 
form; 


^ am  ^ mm  mm  mm  tm  mammmmmmmmmmmmmmmfmmamammmmm-^mmammm’m>mm»mmmt 

(3)  Ho  Bo  Salzer  and  R,  Biicker,  "Tables  of  the  Zeroe  and  Wei^t 
Factors  of  tho  First  Fifteen  laguerre  Polynomials",,  Bull.-. 
Amero  Ifeth  Soc.,  ^55)  1949,  ppol004-1012o 


ts 
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oC, 


abscissa  for  I«guerre  quadra turo* 
wei^t  for  loguerre  quadrature} 
abscissa  for  Hermite  quadrature} 


e weight  for  Hermite  quadrature# 
The  quemtity  Q is  defined  in  Appendix  1, 


It  will  be  observed  that  with  this  formulation,  it  is 

necessary  to  calculate  each  Bessel  function  for  Mi  values  of  its 

KL 

argument.  If  we  assume,  os  above,  that  the  constants  i;  ▼anlsh 

kl 

unless  |k-L|c  3*  then,  in  each  case,  seven  values  of 

k and  of  1 must  be  considered.  With  this  range  of  k £uid  1,  and  with 
B « M = 6,  a total  of  about  3,000  Bessel  functions  Is  required  for 
each  transition.  This  would  be  prohibitive  without  fairly  aimple 
approximations  to  the  Bessel  functions  in  the  range  of  interest# 


For  the  hydrogen-oxygen  reactions,  the  Bessel  functions 
^ p (x)  occur  with  paracaters  rougiily  in  the  range 

p^  12,  V2<~<2# 

t 

Asyn^stotic  formulas  for  Bessel  ftmctions  of  large  order  and 
argument  ere  well-known ^4),  both  for  the  case  x<  p and  x>p#  The 
accuracy  of  the  formulas  decreases  as  x *=>  p,  however,  and  rather 


(4)  G,  B#  Watson,  "Bessel  Functions",  Cambridge,  1945,  Chap,  VIII 


than  add  higher  order  approxlssationa  we  have  found  It  expedient  to 
develop  new  formulae  for  the  region  z ‘^p.  By  eon^ring  the  errora 
in  the  various  approxinationa,  we  have  found  that  the  standard  asjonp- 
totic  formulas  could  be  used  In  the  reglone  x/p<  0o85>  x/p>lol5, 
and  the  new  formula  In  the  region  0,85<^  x/p<1.15,  with  a maximum 
error  of  about  five  percents  The  largest  errors  occur  in  narrwj 
bands  about  the  values  jc/p  = 1 + 0^15»  and  for  only  the  smaller  values 
of  p (p»12),  with  the  er^nrs  decreasing  sharply  away  from  thes^^ 
points..  It  is  felt  that  errors  in  the  integrations  resulting  from 
these  approximations  should  not  exceed  about  two  percent  ^ The  formu- 
las used  in  the  confutation  of  the  Bessel f unctions  are  given  in 
Appendix  II, 

The  properties  of  the  quantity  assumed  above  enable 

US  to  make  a si'dsetantial  reduction  in  the  number  of  Bessel  functions 


that  need  be  computed^  The  orders  of  the  Bessel  functions  have  the 
foxmi  + k(jp»-l)  ^era  C is  a large  constant  (C  > I50,  say)  and 
since  vs  na\st  have  jk-L}A3»  there  is  only  a small  percentage  varia> 
tlon  in  the  oxters  within  a given  transition » Ind«»ed,  as  k ranges 


from  L to  L + 3g  the  order  of  the  Bessel  function  will  Increase 
by  less  than  ten  pereent.  This  suggests  that  for  jlt-lji 
can  be  obtained  by  means  of  an  appropriate  expansion  about  1 

We  have  developed  approximate  formulas  for  calculating  (x) 

in  terms  of  ^p(x)  idiich,  in  the  range  of  Interest  are  in  error  by 
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at  Bost  about  five  percent  for  p with  the  enror  decreasing 
sharply  with  increasing  p.  It  is  felt  that  this  approximation  should 
result  in  errors  not  exceeding  about  two  percent  except  possibly 
for  the  lowest  rotational  states  where  the  error  might  be  as  Mgb 
as  five  percents  The  approximations  are  given  in  Appendix  II. 

It  will  be  observed  that,  with  this  ainplificatlon,  only 
one>'Seventh  of  the  Bessel  functions  need  be  coo^jubed  directly  from 
the  asymptotic  formulas  given  above.  Since  a large  portion  of  the 
computing  tins  is  devoted  to  the  calculation  of  Bessel  functions;; 
this  results  ir  a major  reduction  In  machine  time  for  the  problem. 
Indeed,  without  some  approxiiaation  of  this  sort,  the  conqnittng  time 
for  the  procedure  outlined  here  would  be  prohibitive.. 

Remarks  cc  Coding 


The  work  of  coking  the  present  calculation  for  SE&C 
begins  with  the  fonmilaLlou  presented  here.  The  major  diffieulxies 
encountered  in  the  coding  have  been  chiefly  concerned  with  storage 
limitations  of  the  machine.  In  order  not  to  exceed  the  machine 
capacity,  it  has  been  necessazy  to  devise  several  flexible  oalciilat- 
ing  routines  that  could  be  used  for  more  than  one  purpose  , Th\is, 
for  example,  a single  routine  has  been  set  up  for  calculating 


in  both  regions  x/p^  O.S5  and  x/p  ^ 1,15=> 
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The  developaeot  of  these  procedures  has  been  carried  on 
Jointly  by  representatives  of  hydrocarbon  Research,  Inc„  and  the 
Conputation  laboratory.  At  the  present  tln»,  all  of  the  major 
planning  probleas  have  been  resolved,  and  the  remainder  of  the 
vorii:  at  the  Computation  Laboratory  will  be  ooncemed  with  putting 
the  various  segments  ihto  a contiguous  unit  £uid  with  final  checking 
of  the  code  on  SEACo 


Tho  phsTsical  constants  eniploysd  is  tho  calculation  of  the 
treuisltion  probabilities  an  defined  as  follows!! 
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m Boltzmann  ^8  constant^* 

>=  oqulllbrluB  distance  between  atoiDS  of  initial  molecule 
e equilibrium  distance  between  atoms  of  final  Boleoule.> 

" fundamental  vibration  frequency  of  the  initial  molecule^ 
= fundamental  vibration  frequency  of  the  final  molecule. 


J ^ = moment  of  inertia  of  initial  molecule., 

» aaas  of  final  atom  , 

» mass  of  initial  at on  > 

B^  Bass  of  exchanging  atoB. 
t «=  tiae 

^ = region  of  norBaIi2;ation  of  the  wave  functions  ^ 


The  nuaerical  values  of  the  constants  for  the  hydrogeo'  - 


oxygen  reactions  are  given  below? 


©2  > a--- »0H  -f  0 * 0 — ^ OH  + E 
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There  Is  considerable  uncertainty  in  the  form  to  be 
ployed  for  the  chemical  per1;urbetlon  potential  7^( ^ » We  ha-?® 
found  it  conTenient  to  fit  this  potential  in  the  form 


Y -- 


— I X i / 


! t/- 1^1 

Here,  3 and  ^ denote  the  vectors  joining  th&  atoms  of  the  Initial 
and  final  molecules,  respectively:  In  this  form,  it  will  be  noted 

that,  for  large  R,  the  buUc  of  the  contribution  comes  in  the  neigh- 
borhood of  X = - 1,  corresponding  to  a linear  configuration  of  the 
three  atoms  <.  The  extent  to  which  the  contribution  Is  limited  to 
this  region  depends,  of  course,  on  the  magnitude  of  R,  and  It  would 
be  valuable  to  explore  the  effect  of  changes  in  R upon  the  transl'- 
tlon  probabilities.  Preliminary  considerations  in  the  case  of  the 
l^rogen~oxygen  reactions  suggest  that  we  take  R = 5,  in  which  case 
the  expansion  of ^ in  a Isgendre  series  leads  to  the  following  values 
of  the  coefficients s 

- o85,  ii  e=  ♦ ,52,  h « - .28, 


&^5 


If  we  take  R •=  10,  ue  fixid 

“ -»90,  J2  = «.,73  = -.53 

^ ^ 

Ko  attenqpt  has  been  toade,  at  this  st^ge,  to  evaluate  the  coxiatant  He  j 
Ita  value  will  affect  only  the  abeolut.e  react-loE.  rate  and  not  the 
relative  values  of  the  transition  probabilities. 


Annendl-r  II  . for  BeggeX  Functions 

The  fonaalas  for  oonq^titlng  the  Beesel  functions  are  as 


follows^  let  i - I ' 
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and  let 
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For  7 < ^B5, 
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For  7 > lol5  j 
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The  foimuiab  for  calculating  the  effect  of  small  changes 


in  the  ordei*s  of  the  Bessel  functions  are  as  follows ; 


For  y<.85^ 


For  y > 1,15 
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